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Abstract 

We analyze recent results of SU(3) lattice QCD calculations with a phenomeno- 
logical parametrization for the quark-gluon plasma equation of state based on a 
quasi-particle picture with massive quarks and gluons. At high temperature we ob- 
tain a good fit to the lattice data using perturbative thermal quark and gluon masses 
from an improved HTL scheme. At temperatures close to the confinement phase 
transition the fitted masses increase above the perturbative value, and a non-zero 
(but small) bag constant is required to fit the lattice data. 

1 INTRODUCTION 

Strong interactions are described by SU(3) Yang- Mills field theory, and the fundamental 
degrees of freedom of Quantum Chromodynamics (QCD) are gluons and quarks. At 
high temperature this theory is weakly coupled, allowing for the use of perturbative 
methods. The leading order contributions to the characteristic collective excitations and 
the equation of state of a quark-gluon plasma were already determined many years ago 
^ ^ . A gauge invariant approach to non-leading corrections was derived more recently 
in form of the Hard Thermal Loop (HTL) approximation 0, |, |^ and its improved versions 

The theory of QCD predicts the appearance of a phase transition between the quark- 
gluon dominated high energy region and the hadronic state in the low energy region. 
The two states are characterized by a dramatic difference in the number of degrees of 
freedom. Perturbative QCD can describe successfully only the asymptotic state at very 
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high momenta, but it fails close to this phase transition where non-perturbative effects 
become dominant. Where the perturbative region begins is still a matter of debate 0]. The 
numerical method of lattice QCD can describe both sides of this phase transition. Recent 
developments of this field have yielded dramatic improvements both in the extrapolation 
to the continuum limit and for the inclusion of dynamical fermions. 

The existence of this phase transition generated a large experimental effort to create 
it in the laboratory and investigate it in detail through heavy ion collisions. However, 
to verify the appearance of a deconfined state in such experiments we need accurate 
knowledge about its subsequent hadronization. There is a general interest to create phe- 
nomenological models for this phase transition which agree with the perturbative results 
at high energy and with the lattice QCD data at low energies. Such models must specify 
the basic degrees of freedom in the plasma state which will participate in the formation of 
hadrons during the confinement phase transition. The experience from phase transitions 
in solid state physics and other fields suggests the introduction of quasi-particles with ef- 
fective masses generated through the interactions among the basic constituents. If a large 
part of the interaction can be included into the effective masses, then a quasi-particle 
description of the interacting matter can be generated in which the quasi-particles move 
freely or interact only weakly |[T^. Such a model would be a preferable starting point for 



phenomenological hadronization models and for a description of the strongly interacting 
matter near the phase transition. 

Our goal is to identify the appropriate degrees of freedom through an analysis of 
lattice QCD data. Similar attempts started as soon as good lattice QCD results on the 
equation of state of strongly interacting matter became available, and they developed 
in parallel with the improvements of the lattice data. It was observed already quite 
some time ago that the formula for the free energy density of a massive ideal gas gives a 



quite satisfactory description of the numerical lattice QCD simulations [jTT|. In Ref. |12 
temperature dependent screening masses and coupling constants were extracted from 
SU(2) lattice data [|13[ and compared with perturbative results. In Ref. [|14| it was shown 



that the obtained thermal mass can be indeed used as an effective mass in the equation of 
state, comparing favorably with SU(2) lattice QCD data. Early (very poor) lattice data 
on the pressure and energy density of pure SU(3) gauge theory in the temperature region 



1.2 < T/Tc < 2.4 were interpreted with a constant gluon mass, Mg ^ 500 MeV, 
and a constant bag constant B^^^ «i 200 MeV [0. Newer SU(3) lattice data of better 



quality |]T^ made a reanalysis of the gluonic equation of state possible which yielded a new 
expression for the temperature dependent coupling constant in the thermal gluon mass 
P^ . A quasi-particle based, thermodynamically consistent analysis of the SU(2) |T^ and 
SU(3) [1^ lattice data was performed in Ref. |T0| whose authors investigated in detail 



the temperature dependence of the thermal gluon mass and of the bag constant. New 
SU(3) lattice data with a complete continuum extrapolation appeared in Refs. |]T9|, ^ 



and again a phenomenological analysis proved the applicability of an equation of state 
with massive gluons ||21|. Furthermore it was shown |^ that at T > 2.5 Tc the Debye 
screening mass extracted from lattice QCD correlation functions and the thermal gluon 
mass fitted to the lattice QCD equation of state are consistent with each other, such 
that at high temperature the perturbative QCD with effective massive degrees of freedom 
provides a good description of the lattice QCD results PT| . In Ref. [PTI] the lattice data fit 
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yielded d = 17.2 gluon degrees of freedom, supporting the presence of massive transverse 
modes only. A direct investigation of the screening mass in SU(3) lattice calculations was 



performed in Ref. [E^ on both sides of the phase transition. 



A somewhat different analysis of pure SU(3) lattice data was presented in Ref. |Tj 
It assumes the existence of massless gluons above a certain minimal momentum ("cut- 
off model" ) together with glueball-like non-perturbative massive excitations. This model 
could reproduce the SU(3) lattice data quite well fl^, p3| . 

All of the above analyses were performed for the case of pure SU(2) and SU(3) 
gauge theories, because sufficiently high-quality lattice QCD data existed only in these 
cases (for SU(3) the latest results can be found in Refs. |19|, ^). However, new lattice 



data including dynamical fermions appeared recently, with Nf = 2 [^, ^ and Nf = 4 



26[] quark flavors. Extended investigations are under evaluation for the cases Nf = 1,3,6 
and Nf >7 (see Ref. |2^), but in their present state they can not yet be used to extract 



temperature dependent thermal masses. These new lattice results can now be used for 
the extraction of a phenomenological equation of state for the quark-gluon plasma (QGP) 
containing massive gluons and quarks, which to our knowledge has not been published 
elsewhere. Similar work is under way in the group of A. Peshier and B. Kampfer p8[ . 



The understanding of the dynamical generation of effective quark and gluon masses 
could be very important in many research fields. Dynamical mass generation is an essential 



ingredient for the solution of the infrared catastrophe in hot gauge theories |29] and the 



formulation of effective field theoretical approaches to thermal QCD |^0[. It is at the 



heart of recent new approximation schemes such as the "screened perturbation theory'' 



3T|, |32|. Screening masses have also successfully been used to regulate the low-p^ behavior 
in the parton cascade approach to ultrarelativistic heavy ion collisions |Q. Massive 
quarks are the basic degrees of freedom in hadronization models: the phenomenological 
model ALCOR is based on the coalescence of massive quarks and antiquarks into hadrons 



3^ ; transport descriptions of the hadronization based on the Nambu-Jona-Lasinio model 
contain massive quarks as well [^. Massive quarks and gluons obtained from HTL 
approximation were already applied to estimate charm production at RHIC and LHC 
energies [^. 

In Sec. ^ we summarize the present knowledge about the thermal and screening 
masses in perturbative thermal QCD. In Sec. || we overview the equation of state of 
quasi-particles, its most important thermodynamical properties and discuss the criterium 
of thermodynamical consistency. In Sec. ^ we analyze pure SU(3) lattice QCD results; we 
determine the relevant number of gluonic degrees of freedom, the temperature dependent 
gluon mass and coupling constant, the temperature dependent bag constant and investi- 
gate the screening mass. In Sec. ^ we repeat this analysis on lattice data with Nf = 2 
and Nf = 4 dynamical fermion species. In Sec. § we generate the equation of state of 
the QGP in general and investigate the speed of sound in this system. In Sec. |^ we in- 
clude hadronic matter into our investigation and describe phenomenologically the phase 
transition between QGP and hadronic matter on the basis of our new equation of state 
obtained from lattice QCD data. In Sec. || we discuss our results. 
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2 THERMAL AND SCREENING MASSES IN 
PERTURBATIVE QCD 



2.1 Gluons 



There are essentially two ways to define an effective dynamical gluon mass: via the pole 
of the effective gluon propagator, or via the long-range behavior of the potential between 
two heavy color sources. In both cases one has to address the question of gauge invariance 
of the result for the mass, since the defining objects are not themselves gauge invariant. 

On a perturbative level, it was noticed by Klimov and Weldon that the 
leading term in a high-temperature expansion of the 1-loop gluon polarization tensor was 
gauge invariant. It was later shown by Heinz |37] that the same result for the gluon 



polarization tensor could be obtained from classical color kinetic theory in the linear 
response approximation. From this expression for the polarization tensor Weldon [Q] 
derived the following dispersion relation for transverse gluons with momenta uj,k ^ gT: 

cu' = e + Ml^{T), (1) 

MlJT) = 9^ + ^) . (2) 

One might object that the high temperature limit in which the polarization tensor was 
obtained is inconsistent with the limit uj,k ^ gT under which Eq. (0) was derived. 
However, Rebhan et al. showed recently 0, |[ that the vahdity of Eq. (|I]) does not 
depend on the high temperature limit taken in Ref. 0; indeed, near the light cone the 
transverse 1-loop polarization function is given by the gauge invariant result 

U,{u;' = e) = Ml^iT) (3) 

independent of the magnitude of u = k relative to T. Since for large momenta k ^ gT 
the gluon dispersion relation moves arbitrarily close to the light cone, Mg{T) can be 
interpreted as the thermal gluon mass in the high momentum limit. Its inclusion as 
a gluon mass term in higher orders of perturbation theory removes certain classes of 
collinear singularities near the light cone, and within the context of Hard Thermal Loop 
resummation ^ ^ it leads to the so-called "improved HTL resummation scheme" 

The "Debye screening mass" M^, on the other hand, is related to the behavior 
of gluonic excitations at small momenta. It can be defined either through the static 
limit of the gluon polarization tensor 'n.^y{ijj = 0,k) for — > or via the behavior of 
the potential between two heavy color charges at large distances, V{r) ~ exp{—Mor)/r 
0. In the first case a gauge invariant result can be obtained from the leading term 
of the high temperature expansion for the 1-loop gluon polarization operator 0, ^ or, 
equivalently, from the 1-loop result within the HTL resummation scheme [§, ^ § for the 
gluon polarization operator at momenta u,k < gT. In the second case, care must be 
taken to make the external sources gauge invariant 0. 

In the HTL approximation one obtains from the behavior of 
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Ul{uj = 0,k-^0) the Debye mass 0, |, |8[ 

Ml=2Ml^{T)=g^T'[^ + ^) . (4) 

This result is related to the 1-loop plasmon frequency, Up, obtained from the gluon dis- 
persion relation (pole of the propagator) for /c — > 0, by M|, = 3ujp. 



2.2 Quarks 

Similarly to the gluon propagator, one can study the quark propagator in the 1-loop 
approximation. In the high temperature or low- momentum limit iJ,k <^ T p|, |^ or, 
equivalently, in the HTL approximation ^ one obtains 

Mj(T) = — = . (5) 

This can be interpreted as the effective quark mass for soft quarks with momenta uj,k < 
gT. For high momenta u!,k ^ gT the fermion dispersion relation again approaches the 
light cone, and one can make use of the gauge invariant light cone limit of the 1-loop 
fermion self energy |Q: 

J:f{u;' = e) = M|^(T), (6) 
Ml^iT) = 2MJ(T) = ^. (7) 

Again this is an exact 1-loop result, independent of the value of u = k, and Mf^^o can be 
interpreted as the thermal quark mass for high momentum quarks, k 3> gT. 



2.3 Temperature dependent coupling constant 

The thermal gluon and quark masses obtained from perturbative QCD are displayed in 
Eqs. d^) and (|^). All masses depend on the temperature T and the strong coupling 
constant g. At higher orders of the loop expansion, the coupling g begins to run as a 
function of T, giving rise to a temperature dependent effective coupling constant g{T). In 
SU{Nc) gauge theory at T = 0, in the presence of Nf quark flavors the 1-loop expression 
for the running coupling constant as a function of the momentum transfer Q is 

^'^^'^ ^ (llAr,-2Ar^)ln(-gVA^) " 



Here A^r is the cut-off parameter. It was shown in Ref. that in a thermal system 
it makes sense to introduce a temperature dependent function g^{T) by the following 
parametrization: 

9\T) = ,,,,, o..M^^..^^ .^ • (9) 



{llN,-2Nf)lnF^{T,T„A) ' 
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We expect F(T, T^, A) to be linear in the region where perturbation theory is vahd |39 . 
In other words, writing 

F(r,T„A) = K(T/T,)-^--^ (10) 

we expect K{T /T^ to be a constant in the perturbative region, i.e. at high T. At low T 
a possible T-dependence of K reflects non-perturbative corrections. 

In SU(3) at different Nj we have different critical temperatures Tc and different 



normalization parameters T^/ Aj^ |^ 



T 

S\J(3),Nf = — > T^ = 260 MeV; — ^ = 1.03 ± 0.19 ; (11) 

SU(3), Nf = A — ^ Te = 170 MeV; = 1.05 . (12) 

For Nf = 2 the critical temperature in SU(3) was determined as = 140 MeV [^, p5|] , 
but the normalization factor Tc/Ajj^ is not known by the authors. From Eqs. ([ll|) and 
( p!2D we will estimate this value by interpolation and we will use Tc/Ajjg = 1.03. 



We will try to fit the lattice data with a non-interacting gas of massive quarks 
and gluons, with masses given by the perturbative expressions (|]) and (0), but with g"^ 
replaced by a phenomenological running coupling g^{T). Using Eqs. (^) and ( p!OD we can 
then extract a function K{T /T^. In a very early work [|12[ in pure SU(3) gauge theory 
this function K{T /T^ was determined from a similar fit to the heavy quark-antiquark 



potential; the fit result was a constant, K{T/Tc) = 19.0. However, the authors of Ref. [12 
used Tc/ Ams = 1-78 ± 0.03 which differs strongly from the now accepted value. With our 
normalization (pJ]), K{T/Tc) = 19.0 would correspond to K{T/Tc) = 33.8. 



Another result was obtained recently from a numerical fit of the equation of state 
in pure SU(3) gauge theory [pT| : 



F(T, T„ A) =4.17^-2.96 . (13) 



Obviously, this parametrization differs from Eq. (p^O]). 

Our expectation is to obtain a constant K{T/Tc) function in that temperature re- 
gion where perturbative QCD and the quasi-particle picture of the lattice QCD results 
overlap. A strong temperature dependence of K{T /T^ implies non-perturbative effects. 
We summarize our results on the function K[T /T^) for different values of Nj in Sects. ^ 
and ^. 



3 EQUATION OF STATE WITH QUASI- PARTICLES 

In this Section we introduce the equation of state of a pure gluon plasma and a quark-gluon 
plasma consisting of massive quasi-particles. This equation of state (resp. its parameters) 
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will be determined from the SU(3) lattice data. The key point is to consider temperature 
dependent effective masses, Mj(T), which are dynamically generated. We can introduce 
the dispersion relation for quasi-particles: 

u;^{k) = e + Mf{T) . (14) 

Here k is the momentum of the quasi-particle and uJi{k) is its energy. A pure gluon 
plasma contains only gluons {i = g, Nj = 0), a QGP consists of quarks and antiquarks 
also {i = g,q,q, Nf > 0). We can define the Bose or Fermi distribution functions fi{k) 
for the quasi-particles as 

Mk) = [exp{^k^ + Ml{T)/T} ± 1]-' = [exp{uj,/T) ± 1]-^ . (15) 

In addition to the effective masses Mi{T) we also need the effective number of degrees 
of freedom of the quasi-particles to determine pressure and energy density. The correct 
value for the number of effective gluonic degrees of freedom is not obvious, since a massless 
vector boson has two helicity states, but massive vector bosons have three spin states. 
One possibility is to introduce a temperature dependent effective degeneracy factor D(T) 
for the gluons. We will determine the value of -D(T) in a pure gluon plasma {Nf = 0), 
where quark degrees of freedom do not interfere. Later we can assume the same value for 
D{T) in a QGP. 

Another possibility is to fix the number of gluonic degrees of freedom (e.g. at its 
perturbatively expected value Dg = 16) as well as the number of quark and antiquark 
degrees of freedom. In this case we have the freedom to introduce an effective interaction 
term for those contributions of the strong interaction which can not be absorbed by the 
presence of effective masses. We parametrize this term by a temperature dependent "bag 
constant", B{T). This interaction term B{T) can be determined in a pure gluon plasma 
as well as in a quark-gluon plasma. As we will see in the following Section, in the QGP 
phase quarks and anti-quarks will also contribute to B{T), so its physical meaning differs 
from the bag constant parameter of the MIT Bag Model; nevertheless we will use the 
same notation for convenience. 

3.1 The number of effective gluonic degrees of freedom D{T) 

Let us first consider pure SU(3) gauge theory where we wish to find out the number 
of contributing gluonic quasi-particle degrees of freedom, D(T). A strongly interacting 
gluon plasma contains not only the transverse gluon modes, but at low momenta k < gT 
there exist also longitudinal plasma excitations. It was shown in Ref. that at high 
momenta k ^ gT the longitudinal modes disappear, in the sense that the residue of the 
corresponding pole in the propagator becomes exponentially small. Since the equation of 
state is dominated by particles with momenta A; ~ T, we expect at high temperatures T — > 
oo (where g becomes small) the contribution of the longitudinal modes to be negligible. 
On the other hand, at low temperatures (where g can become of order 1 or larger) their 
contribution may be relatively large. Therefore we can not a priori neglect them from the 
analysis of lattice QCD data. 
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We will introduce a temperature dependent effective number of gluon degrees of 
freedom D{T), generalizing the approach of Ref. where D was taken as a constant 
(temperature independent) fit parameter. We expect that in the high temperature limit 
there are only transverse gluonic quasi-particle modes, yielding D[T = oo) = Dg = 16. At 
low temperatures the longitudinal gluons are expected to also contribute, giving D{T) > 

With these assumptions the thermodynamical pressure Pg{T) of pure gluon matter 
(which is nothing but the grand canonical thermodynamical potential) can be written as 

From the thermodynamic identity 

e{T) = T . - P{T) (17) 

(we work at zero net baryon density, fiB = 0) one obtains the following relation for the 
energy density: 



n(T) roo 



rr.P'AT) dD{T) ,^,dMJT) D(T) d^k ,,,, 

Here Wd{T) summarizes the extra terms stemming from the temperature derivative of 
the effective mass and effective number of degrees of freedom. On the other hand, a 
consistent quasi-particle picture demands Wo = ||lOl- Requiring this identity yields a 
differential equation for D{T): 

D{T) = DiT*) + r dT'MgiT')^^^ 3DiT') • (20) 

Here T* is an integration constant. If the functions Mg{T) and D{T) satisfy this selfconsis- 
tency condition, then we have a thermodynamically consistent quasi-particle description. 

Lattice calculations yield numerical values for the functions Pg{T) and Sg{T). These 
two sets of data are sufficient to determine for each value of T the quantities Mg{T) and 
D{T) such that they satisfy Eqs. (|TB|-|T^) with Wn = 0. Since the lattice results were 
created in a thermodynamically consistent way we expect that the resulting functions 
Mg{T) and D{T) will satisfy Eq. (^Of) . We will check this as a test for the consistency 
of our extraction procedure for Mg{T) and D{T). Our numerical results on Mg{T) and 
D{T) will be summarized in Sec. |. 



3.2 The temperature dependent interaction term B{T) 

We can also investigate the high temperature limit from another point of view [|1^, Let 
us fix the value of the gluonic degrees of freedom at every temperature to D(T) = Dg = 16, 
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assuming that only transverse gluons contribute to the thermodynamic quantities, and 
calculate directly any extra contribution to the pressure and the energy density. In this 
way both pure gluon matter and a quark-gluon plasma can be investigated. We will here 
consider a QGP with Nf quark flavors and the usual number of quark degrees of freedom: 
Dg = Dg = 2-3- Nf = QNf. 

The effective quark mass Mq{T) will be related to the effective gluon mass Mg{T) 
through the value of the temperature dependent coupling constant g(T) via the pertur- 
bative expressions in Eqs. (|^) and for the effective masses. This means that we absorb 
non-perturbative features into a common non-perturbative flt function g{T), without 
touching the perturbative form of the Eqs. @) and itself. 

The deviation of So{T) = J2i^i(T) and Po(T) = J2iPiiT) from the ideal gas values 
corresponding to the effective masses Mi{T) and the flxed degeneracy factors Dg, Dg and 
Dq will be parametrized by a temperature dependent function B{T) . At high temperature, 
where we expect a perturbative picture based on free quarks and transverse gluons to be 
valid, B{T) should vanish. At lower temperature, B{T) provides another measure for non- 
perturbative physics which can not be absorbed into effective quark and gluon masses. 
Both quarks and gluons contribute to B(T). 

We introduce B{T) by writing the pressure in the form 

r). POO u'^ 

PiT)= E 7^/ d'k-—Mk)-B{T), (21) 
iJ^^g (27r)3 7o 3ui{k) 

motivated by the MIT Bag Model. Thermodynamic identities yield the following relation 
for the energy density: 

^(T) = E 7^ / d'k u,{k) m + B{T) + WsiT) . (22) 

Here Wb{T) summarizes the extra terms obtained from the temperature derivative of 
B{T) and of the effective mass Mi{T). Similarly as above thermodynamic consistency 
of our quasi-particle model requires Wb{T) = 0. This yields to the following integral 
representation for B{T) W^, 



B{T) = B{T*)- E tAi/ dT'M,{T'f-^^ —Mk). (23) 



i=g,q,q 



(27r 



Again T* is an integration constant. If this expression is satisfled by the functions Mi{T) 
and B{T), we have a thermodynamically correct quasi-particle description containing 
effective gluons and quarks with effective thermal masses. 

Similarly as above, from the numerical values for the functions P{T) and e{T) 
obtained from lattice QCD and assuming Wb{T) = 0, we can determine the appropriate 
values for Mi(T) (resp. g(T)) and B(T) through Eqs. (^) and (|22|). As before we expect 
that the obtained functions Mi{T) and B{T) will satisfy Eq. (pS)). Once again this is a 
good numerical consistency check for our extraction of the functions Mj(T), B{T). 
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4 ANALYSIS OF LATTICE DATA FOR PURE 
SU(3) GAUGE THEORY 



We analyze the latest lattice data in the continuum limit for pure SU(3) gauge theory 
|T9| , |20[| . The calculations of Ref. ||19[ show a ±2% error on the data for -P(T), e{T) and 



s(T). We include this error into our analysis; it will be indicated in Figs. 1, 2 and 3 by 
dotted lines above and below the solid lines corresponding to the results extracted from 
the mean values. 

We begin by determining the effective number of gluonic quasi-particle degrees of 
freedom in a gluon plasma, extracting the functions Mg{T) and D{T) as explained in 



Sec. |3TT| . Fig. la shows the temperature dependent effective gluon mass, Mg{T), while 
in Fig. lb we display the effective number of gluonic degrees of freedom, D{T). At 
high temperature, T > 2.3 Tc, the number of degrees of freedom is rather constant, 
D{T) = 14.5 ± 1. Please note that the ±2% error of the lattice data results in an 
uncertainty of about 1 gluon degree of freedom. Within this uncertainty the result is 
surprisingly close to the naive expectation that at high temperature only the transverse 
gluonic modes are present. In the low temperature region, T < 2.3 Tc, the function 
D{T) increases, which one may wish to associate with the increasing contribution of the 
longitudinal modes. However, for T < 1.5 Tc, D{T) rises very strongly and exceeds even 
the value of 24 expected for massive gluons with 8 color and 3 helicity states. This 
indicates the presence of strong non-perturbative effects at low temperature not all of 
which can be absorbed by the effective gluon mass. In fact, part of the strong rise in 
D{T) is driven by the strong increase of the effective gluon mass which leads to exponential 
suppression of the thermodynamical contributions to P{T) and s{T). We consider this 
behavior of D{T) as physically unreasonable and believe that a parametrization of these 
non-perturbative effects via an extra interaction term makes more sense. 

To this end we fix the effective number of gluonic modes to Dg = 16 as suggested by 
the behavior of D{T) in Fig. lb and our perturbative prejudice at large T. This means 
that we consider only transverse effective gluons as quasi-particle states and subsume all 
further interaction effects into a non-perturbative "bag constant" B{T). The procedure 
is described in Sec. |^. Fig. 2a shows the temperature dependent transverse gluon mass 
Mg[T) resulting from this approach. Comparing with Fig. la one sees that at large values 
of T the effective mass is now somewhat larger, compensating for the choice of a fixed 
Dg = 16 instead of the best fit value Dg ^ 14.5 in Fig. lb. For smaller temperatures, 
however, the rise of Mg{T) is much weaker than in Fig. la, and Mg{T) also shows a much 
weaker sensitivity to the statistical error of the lattice data on P{T) and e{T). This 
indicates that this type of parametrization provides a more reasonable fit to the lattice 
data than the one above in terms of Mg{T) and D{T). Please note that the effective 
gluon masses shown in Fig. 2a correspond in the region T > 2 Tc to a coupling constant 
g = \/2Mg(T)/T ^ 1.3 — 1.4, in agreement with other extraction procedures: for example. 



the authors of Ref. |41| obtained from the interquark potential at small distances the value 



a ~ 0.15 which translates into g = 1.37. 

Fig. 2b shows the function B{T)/e{T), the interaction term relative to the total 
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energy density. At high temperatures, T > 2.3 T^, the bag constant is very small, only 
about 2 — 3% of the total energy density, although slightly negative. We take the smallness 
of B(T) as confirmation for the validity of a quasi-particle picture at large temperature. 
Near T = 2.3 Tc B{T) changes sign, increasing for smaller values of T, but never ex- 
ceeding a value of 15% of the total energy density even at Tc. Note that this relatively 
small deviation from an ideal quasi-particle picture gave rise to the dramatic and un- 
physical behavior of D{T) (and Mg{T)) in the alternative approach above. Clearly, the 
parametrization via B{T) is more economical and leads to a more selfconsistent picture. 

The dashed line in Fig. 2b shows the result of an integration of Eq. (|23|) , using the ex- 
tracted Mg{T) from Fig. 2a and choosing boundary conditions B{T*) = at T* = 2.3 Tc. 
The agreement with the numerically extracted values for B{T) is nearly perfect indicating 
the thermodynamical consistency of the lattice data and our extraction procedure. 

From the temperature dependent gluon mass, Mg{T), one can determine the tem- 
perature dependent couphng constant, g{T), and extract the function K(T/Tc) given in 
Eq. (0). The result for pure SU(3) gauge theory is shown in Fig. 3. The dotted lines 
again indicate the ±2% systematic error of the lattice data. One can see that at high 
temperature we obtain a constant value, K{T/Tc oo) = 18 ± 8. At low temperature 
K[T/Tc) decreases indicating a larger coupling constant and a larger thermal gluon mass. 
Close to the phase transition the function approaches the value K{T/Tc — > 1) = 1.5; the 
precise value depends, however, on the latent heat at the phase transition point which is 
still under intensive investigation in lattice QCD calculations. 

The obtained curve can be fitted rather well with the following function (dashed 
line in Fig. 3): 

If we start from this approximate function of K{T/Tc) for the coupling constant g{T), 
we obtain an analytically parametrized, approximate thermal gluon mass Mg{T) (see 
Eqs. (p|,|9| JT0|) ). Once we fix the boundary condition B(T*) = at T* = 2.3 Tc as deter- 
mined from the numerical analysis of the lattice data (see Fig. 2b), the bag constant B{T) 
can be uniquely reconstructed by integrating Eq. (|23|) . With Mg{T) and B{T) known, 
P(T) and e(T) are easily evaluated using Eqs. ( pT|p2D (with Wb = 0). This means that 
Eq. ( P^ together with the value for T* provide a complete analytical parametrization of 
the equation of state. In Fig. 4 we compare the original lattice QCD data (diamonds) 
with the thus reconstructed values (dashed and dash-dotted lines). The dash-dotted lines 
neglect the interaction term B{T) and thus represent only the quasi-particle contribution 
to P and e. Clearly the agreement between the model and the lattice data is at most 
qualitative in this case. Including the bag term, however, (dashed lines) we obtain nearly 
perfect reproduction of the lattice data. Fig. 4 demonstrates the usefulness of a thermo- 
dynamically consistent quasi-particle picture. The deviation from the Stefan-Boltzmann 
values indicates that the massless degrees of freedom can not provide a satisfactory de- 
scription, while the quasi-particle model plus bag term leads to qualitative improvements. 
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5 ANALYSIS OF SU(3) LATTICE DATA FOR 

A^/ = 2, 4 



Newest lattice calculations including dynamical fermions can be used to determine the 
equation of state for a realistic quark-gluon plasma including light quarks. The numerical 
lattice QCD results can be found in Refs. |2§ for Nf = 2 and in Ref. for Nf = 4. 



Following these articles we interpolated the expected continuum result in both cases. 



From Ref. |2J] for Nf = 2 we used the data for a current quark mass aruq = 0.0125 



(octagons in Figs. 9 and 10 in Ref. |2^ and circles in Fig. 7 here). The statistical error 
of the lattice data is small, < 1% for the pressure and ^ 2 — 4% for the energy density. 
From Ref. |2^ for Nf = 4 we used the data set for a current quark mass mq/T = 0.2 



and we considered their extrapolation of the energy density to the chiral limit (see □ in 
Fig. 4 of Ref. PB| and in Fig. 7 here). In this case the statistical errors are even smaller, 



but the systematic errors resulting from the extrapolations to the chiral and continuum 
limits are hard to judge. To estimate the possible influence of errors on the lattice data 
for our extraction procedure we will consider a universal ±2% error in both cases; the 
corresponding uncertainty will be indicated by dotted lines in Figs. 5 and 6. 

In the framework of perturbation theory the effects of dynamical quarks on the 
gluon dynamics are small and can be simply included into the thermal gluon mass and 
Debye screening length. No additional gluonic collective modes arise. However, at small 
momenta there are additional fermionic collective modes, the "plasminos" [^. Based on 
our experience in pure gluodynamics they are not expected to contribute to the equation 
of state at high temperature; again the residues of their poles in the quark propagator 
vanish exponentially for moments k ^ gT 0]. We will therefore assume a fixed number of 
degrees of freedom, given by the perturbative values at high T, Dg = 16, Dq = Dg = 6Nf, 
and include all collective interaction effects into an interaction term B{T) as for the purely 
gluonic system above. 

The obtained results are qualitatively very similar to the previous ones for the purely 
gluonic case. Fig. 5a shows the obtained results for the effective gluon mass, Mg/T. Since 

^^^^ = ^^y'l + Fig. 5a can be viewed as giving the temperature dependence of the 
effective coupling constant g{T) which according to Eq. (^ also determine the effective 
quark mass. From Fig. 5 and Eq. (^ we extract g ^ 1.6 in the region T > 2 for 
Nf = 4, somewhat larger than for the purely gluonic case. 

Fig. 5b shows the interaction term relative to the total energy density B{T)/e{T) 
both for Nf = 2 and Nf = 4. Similarly to the result in pure SU(3) gauge theory, the 
interaction term B{T) remains small at T > 2 — 2.3 T^, and does not exceed 15% of the 
energy density in the low temperature region. In the case Nf = 2 the lattice data allow 
us to extract Mg{T) and B(T) only in the temperature region Tc < T < 1.4 Tc, however 
the characteristics of the obtained curves are close to the results obtained for Nf = 4. 

Fig. 6 displays the obtained functions K{T/Tc) (solid lines) together with the ex- 
pected errors (dotted lines) stemming from an assumed typical ±2% uncertainty in the 
lattice data. For all three cases, Nf = 0, 2 and 4, the functions K{T/Tc) are very similar 
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(within the admittedly large error bars). Within the existing uncertainties we can thus 
define a universal function K{T/Tc); differences due to variations of A''^ and Nj can be 
largely absorbed into T^. and A^g^ in Eq. (]T0[). It would be very interesting to see whether 
the existence of such a universal function K{T/Tc) is confirmed by future higher quality 
lattice data. 

Since the data for purely gluonic SU(3) span a much larger range than those includ- 
ing dynamical fermions, we will use the analytical form (p4D extracted from the gluonic 
case also in our quasi-particle picture for the full QGP. For the zero point of B{T) we will 
again choose the value T* = 2.3 determined in the purely gluonic SU(3) case, although 
(modulo large numerical uncertainties) the Nf ^ data seem to favor a somewhat smaller 
value. 

For a QGP with N j = 2 quark fiavours Fig. 7a shows the lattice data (circles 
for pressure and energy density) and the reconstructed values for e{T) and P{T), both 
including the bag constant (dashed lines) and neglecting the interaction term (dash-dotted 
lines). For Nf = 4 quark flavours the analogous quantities are displayed in Fig. 7b. Here 
the squares give directly the lattice data from Ref. for the energy density, while the 
diamonds represent the continuous line for the pressure given in Fig. 2 of Ref. ||2^ . In each 



case inclusion of the interaction term B{T) improves the picture. The agreement is weaker 
than in the pure SU(3) case (see Fig. 4), but the pressure is reproduced rather nicely as 
well as the energy density at high temperature. Note that here we used the function 
K{T/Tc) from Eq. (p^). The discrepancy in the region T Tc may be connected to the 
use of massive current quarks in the actual lattice simulations and associated uncertainties 
in the extrapolation to the chiral limit (for a discussion see Ref. [0). In view of those the 
10 — 15% discrepancy between the reconstructed curves and data is surprisingly small, 
especially if we remember that we used here the parametrization which was obtained from 
the pure SU(3) case. 



6 SPEED OF SOUND IN THE MASSIVE QGP 

In hydrodynamical models of strongly interacting matter the dynamical behavior is mostly 
determined by the value of the speed of sound: 

^ dP/dT 

' de de/dT ■ ^ ' 

Since the energy density s{T) was obtained from Eq. (|17D we can rewrite this as 

= (^ + ^)/^ (26) 
de/dT ' 

which is more directly connected to the lattice data. For a noninteracting {B(T) = 0) 
gas of massless particles the speed of sound is = 1/3. In Ref. |3^, using Eq. (p6D, it 
was found that even for a non- interacting {B(T) = 0) massive quark-gluon gas the speed 
of sound remains very close to = 1/3. (In that work thermal masses from earlier HTL 



13 



calculations and a running coupling constant according to Eqs. (|9,1C) with a constant 
K{T/Tf.) were used.) If there are, however, strong remaining interactions which cannot 
be absorbed in the masses but require a non- vanishing and possibly strongly T-dependent 
B{T), this will generate deviations of cj. from the ideal gas value, especially near Tc |0 . 



A strong drop of near Tc will give rise to a "soft point" in the equation of state at 
which the ability of the matter to generate expansion flow is minimal. 

Since in our analysis the interaction term B{T) remains small and changes smoothly 
with T, we expect only weak modifications of the speed of sound around Tc. Fig. 8 
displays the speed of sound for A^j = 0, 2, 4 in our quasi-particle picture and confirms 
this expectation. Including the bag term B{T) (dashed lines), cj. ~ 0.33 in the high 
temperature limit, and decreases to a minimum value chimin = 0.15 around T^. If we 
neglect the interaction term B(T) and calculate the speed of sound g from the expression 

2 _ dPo _ dPo{T)/dT 

deo deo{T)/dT' ^ ' 

where Pq{T) = P{T) + B{T) and eo{T) = e{T) - B{T), then (see dash-dotted line in 
Fig. 8) the speed of sound shows only a very weak temperature dependence, decreasing 
from c^Q ~ 0.30 at T = 4.5 to c^Qlmm = 0.23 at Tc. If we instead use Eq. ( p6D also in 
this case (as in Ref. |^), where B{T) cancels in the numerator, then the resulting speed 
of sound is very close to the full result (dashed line), because of the small influence of 
B{T) on the energy density in the denominator of Eq. (^61) . Our result agrees with the 
calculations in Ref. in the high temperature region and differs around Tc because of 
the non-linear behavior of K{T/Tc) indicating non-perturbative effects. 

These results on the speed of sound show that a massive quark-gluon plasma will 
expand rapidly and that therefore the phase transition should be a fast process without 
detectable duration effects. Even where the pressure already shows large (> 50%) de- 
viations from the massless ideal gas law, the speed of sound still deviates by less than 
10% from the value l/VS. In this way our analysis of lattice QCD results favors fast 



hadronization models [p4l E3l 44] of the deconfined phase. 



7 GLUON SUPPRESSION AROUND IN SU(3) 
GAUGE THEORY 

In this Section we investigate some phenomenological consequences of the large effective 
quark and gluon masses. Fig. 9 displays the effective masses in GeV for Nf = 0, 2, 4. All 
masses behave similarly, since they are connected to each other by the function K{T/Tc) 
in the effective coupling constant g{T). Just above the critical temperature the effective 
masses drop, but then they remain rather constant in the temperature region 1.2 Tc < 
T < 2.5 Tc. The mass values for the case Nf = 2 are smaller than for Nf = due to the 
smaller critical temperature Tc = 140 MeV instead of 260 MeV. For Nf = 4 they increase 
again due to the larger contribution from quark loops and the somewhat larger Tc = 170 
MeV. The gluon masses are large in the cases Nf = 0, 4, of order Mg = 650 — 750 MeV in 
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the temperature region Tc < T < 3 (solid lines). In the same temperature region the 
quark masses are of order Mg = 300 — 400 MeV (dashed lines). 

Large effective masses yield smaller number density at the same energy density. 
Thus the density of these massive quasi-particles will be smaller than in the massless 
case. Fig. 10 demonstrates this by displaying the density ratio of the massive quasi- 
particles and the massless ones, nj/n° (solid lines for gluons and dotted lines for quarks) 
as a function of temperature. At large temperature the ratio flattens and eventually 
approaches 1, but close to the critical temperature the ratio drops very quickly to values 
much smaller than 1. This indicates that the massive quark-gluon plasma is much more 
dilute than a massless QGP, especially near Tc. 

Now we compare the abundances of quarks and gluons relative to each other. 
Eqs. (|^) and (^ imply that in our model the gluon effective mass is always larger than 
the quark effective mass by the ratio 




(2^ 



which even increases with the number of quark flavor Nf. Each gluon degree of freedom 
will thus be suppressed compared to each quark degree of freedom. After multiplying the 
corresponding degeneracy factors we can compute the density ratio for gluons 

= ^""^ , (29) 
ng + ng + rig 

which is shown in Fig. 11 as a function of T. The horizontal lines show the ratio 
for the massless QGP. The deviation caused by the effective mass is clearly seen, and a 
considerable suppression of the relative gluon number occurs near the phase transition. 
Near T^, for Nf = 2 we have Rg = 0.25 which means that only 25% of the particles 
are gluons while the remaining 75% are quarks and antiquarks. For Nf = 4 the gluon 
suppression is even stronger and near Tc we find only 10% gluons and 90% quarks and 
ant i- quarks. 

These results show that the massive QGP is a dilute system dominated by massive 
quarks and antiquarks. Gluons are much heavier than quarks and their number densities 
are suppressed. These properties indicate that lattice QCD results support the formation 
of a quark-antiquark dominated plasma state just above the critical temperature, which 
contains mostly massive quarks and antiquarks. This system hadronizes quickly. Fast 
hadronization is also favored by the fact that the effective quark masses are already close 
to the masses in the "constituent quark" picture of hadrons. Clusterization of these 
massive quarks leads directly to multiquark states with similar masses as those of the 
finally observed hadrons, allowing for immediate hadron formation without having to 
wait for further energy transfer. 
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8 CONCLUSIONS 

Lattice data in pure SU(3) gauge theory and new lattice results with Nf = 2 and Nf = 4 
dynamical fermions suggest that beyond the phase transition temperature strong non- 
perturbative effects dominate the deconfined state. In order to understand this phe- 
nomenon and to facilitate the application of lattice results in dynamical descriptions of 
the evolution of the deconfined phase and its hadronization we analyzed the lattice data 
in a phenomenological model assuming the appearance of massive quasi-particles, namely 
massive quarks and massive transverse gluons. We determined a temperature dependent 
effective coupling constant g{T), which determines the effective masses, and the inter- 
action term B{T), which summarizes contributions not to be included into the effective 
masses. We found that such a quasi-particle picture works very well. The interaction 
term B{T) remains small and the extracted dynamical mass is consistent with perturba- 
tive QCD results at high temperature while it includes non-perturbative contributions at 
low temperature. In all cases, Nf = 0,2,4, the obtained characteristics are very similar. 
If we determine g{T) in the purely gluonic case and extend it in an appropriate way 
for Nf > 0, we can reproduce approximately the Nf = 2,4 lattice data. This fact sug- 
gest the existence of a universal description of non-perturbative effects in the fermionless 
and fermionic cases. Our analysis leads to a thermodynamically consistent quasi-particle 
model for the QGP equation of state which is parametrized via an effective coupling con- 
stant g{T) through Eqs. (P| JT(]|J2^ ) and the zero point T* of the interaction term B{T). 
Further lattice data with improved quality over a wider temperature region are needed to 
fix g{T) and T* more accurately. 

The successful reproduction of lattice data by this quasi-particle model indicates the 
validity of such an approximation. The consequences of the appearance of these quasi- 
particles are very interesting: near Tc lattice results indicate the existence of a QCD 
phase which is dominated by massive quarks and anti-quarks, while the even heavier 
gluons are suppressed. This quark- ant iquark plasma can be characterized by a relatively 
large speed of sound which indicates fast dynamical evolution and the lack of long time 
delays during hadronization. Massive quarks and anti-quarks can form clusters, 2- and 
3-body objects, which can be easily associated with the hadronic mass spectra: massive 
quarks {Mq = 300 — 350 MeV) above the critical temperature turn into similarly massive 
"constituent quarks" inside hadrons below in a smooth but rapid hadronization process. 
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Figure Caption 



Fig. 1: The temperature dependent effective gluon mass, Mg{T), in units of temperature 
T (a) and the effective number of gluonic quasi-particle degrees of freedom, D[T) 
(Jo) in pure SU(3) gauge theory (A^/=0) as a function of T/Tc. The horizontal line 
in Fig. lb indicates Dg = 16. 

Fig. 2: The temperature dependent effective gluon mass, Mg{T), in units of temperature 
T (a) and the interaction term B{T) normalized by the energy density b{T) (b) for 
pure SU(3) gauge theory (A^/=0) as a function of T/Tc. For details see text. 

Fig. 3: The function K{T/T^ for the effective coupling constant g{T) in pure SU(3) 
gauge theory. The solid line shows the mean values, the dotted lines show the 
influence of a ±2% systematic error of the lattice data on the function K{T/Tc). 
The fitted function Eq. (0) is denoted by the dashed line, and the horizontal line 
at = 18 indicates its asymptotic value. 

Fig. 4: Lattice data for pure SU(3) gauge theory on e/T^ and 3P/T^ (diamonds) ||T9(| , 
their reconstruction with a free quasi-particle picture with B{T) = (dash-dotted 
lines), and including the interaction term B{T) (dashed lines). The horizontal line 
at 5.24 indicates the Stefan-Boltzmann value for massless non-interacting particles. 

Fig. 5: The temperature dependent effective gluon mass, Mg{T), in units of temperature 
T (a) and the interaction term B{T) normalized by the energy density e{T) (b) in 
SU(3) gauge theory for Nf=2 and Nf=A as a function of T/Tc. 

Fig. 6: The function K{T/Tc) from the effective coupling constant g{T) in SU(3) gauge 
theory Nf = 0,2,4 (solid lines) and the influence of ±2% systematic error of lattice 
data on the function K{T/T^. 

Fig. 7: SU(3) lattice data (symbols) and our model reconstruction (lines): Nf = 2 |^ 
(a) and Nf = 4 |2^ (b). The reconstructed quantities are denoted by dash-dotted 
lines in a simple quasi-particle picture and dashed lines if the interaction term B(T) 
is included. Horizontal lines indicate the Stefan-Boltzmann values for massless non- 
interacting particles, 12.17 and 19.08, respectively. 

Fig. 8: The speed of sound c^, calculated for Nf = 0,2,4. The dash-dotted lines denote 
the results c^g from the free quasi-particle model {B{T) = 0) while the dashed lines 
with asymptotic value = 0.33 indicate the results including the interaction term 
B{T). Note that the differences between Nf = 0, 2 and 4 are hardly visible. 

Fig. 9: The effective masses of quarks (dashed curves) and gluons (solid curves) in SU(3) 
gauge theory with Nf =0, 2 and 4 quark flavour. 

Fig. 10: The density ratios between massive (nj) and massless (ra?) quarks and gluons 
in SU(3) gauge theory for Nf =0, 2 and 4. 

Fig. 11: The relative gluon density Rg = Ug/ (jig + + rig) for massive, for massless 
particles in Nf = 2 (solid lines) and in Nf = 4 (dashed lines) QGP. 
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